Pederiamos admitir outra situação: 


Que na superfície em movimento P (r, t) apenas 
exiba velocidade tangencial 
v, 


Nestas condições as relações (5b) ainda são 
válidas e como 


Vs dis Ú (57) 


reduzem-se às condições fronteiras para os corpos 
em repouso: 


nX[H] =], (58a) 
nX[E] —- O (58 b) 
n. [D] = p, (58 c) 
n.[B]=0 (58 d) 
n.[3,]=0 (58 e) 
n.[JM] +v,.J,+ == 6 (58 £) 
+ 
* ” 


Uma vez que os resultados obtidos e constantes 
de (56) se suportaram nas leis de transformação 
das densidades de carga e de corrente superficiais 
para um deslocamento segundo a velocidade v nor- 
mal à superficie, a generalização para v £ Vi. exige 
alguma reflexão. 

À fim de esclarecer este assunto consideremos 
S (x.y,z, t) e nesse referencial um plano paralelo ao 
eixo dos zz animado de movimento de translação 
de velocidade 


ve, (59 a) 
Seja n a normal a esse plano no ponto P (x, 
YZ t). 
Seja 
«= <+(n,e,) (59 b) 


pelo que 
n = cosa e, + seng e, (59 c) 


Trabalharemos com a Transformação Especial 
de Lorentz pelo que se identificarão: 


e, =8. e, =e, e“=e, (59 d) 
Admitamos sobre esse plano uma densidade de 
carga superficial p, e uma densidade de corrente 
superficial J,, resultante de uma integração de J e p 
numa espessura ) segundo n. 


42 


De acordo com a atitude tomada na secção 4, 
teremos que no instante t=0 e no referencial 
S' (x, v, 2',t') se observa uma normal 


n' = cosg' e. + seng” e, (59e) 


no ponto P (x',w'z', t') do plano transformado, 

As distribuições J' e p' conducentes a J' e pf, 
resultarão de uma integração na espessura N se- 
gundo n', sendo fácil concluir que 


pa sen q Cos q 
—— — — = y (59 £) 
A sen q COS m 
Assim e atendendo a (8): 
& sm = y [Sm Pv] RPA (60 a) 
Do == Ty NUA (60 b) 
à = Dia NYN (60 c) 
p= Em VT] N/A (60d) 


Posto isto partamos das condições fronteiras 
para corpos em repouso: 


n'X [HH] = Fº (681 a) 
nº. [D'] = pr (61 b) 


Esta condições desdobram-se em: 


cosm [Hº] — seno” [H'] =J'.. (62 a) 

cos [HE] = (62b) 

seng' [H' ]=J'.. (62c) 
acrescidas de 

nt Jo =J. cosq' + é sena! = O (62 d) 
o que aliás está implícito em (62b) e (62c) 

Teremos ainda, de (61 b) 

cosa! [D',] + seng' [DJ = 4º, (646) 


Atendendo a (6), (7), (60) e (59f) teremos: 


| cosa! 
y cosa! | [H,]+v [D, | sena) = Ss o 463 a) 


Ge 


(63 b) 


Y cosa” | LEA] —y (D,1| o ES é 
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sen q 
Sen a 


(63 c) 


y sen q [LH] [D,] |= Y [JPY] 


cos a [D,] + y sena! | [D,|—-— [E] |= 


=y[ p—— M e (6) 
c sen É 

o que simplificado conduz a 

Cos q [ [H,] + v [D,] |— sena [HJ] =J, (64a) 
cosa |— [H,] +v[D,]| =J, (64 b) 
sen x [ [H]—* [D,] = Jp (64 €) 


1 | f v 
—a cosa [D,] + sena | [D,]—-—— [H,] |= 


vI (64 e) 


Salientando-se que (62d) conduz a 
COS q E Pv] + sen q Ea = 0 (64 d) 


Convirá ainda escrever a equação (Gde) sob a 
forma equivalente: 


Y 


cosa [D,] + sena [D,] ——— cosa [D/] — 
af v nm 1 
— Bem q e [H,] DV (64 £) 
Das equações (64) e designando 
cosa v (65 a) 
por 
(v.n) (65 b) 


em que, em (65b), v tem carácter vectorial, resulta: 


cosa [H ]—-sena [H/] +(v.n) [D]=J,  (66a) 
— cosa [H,] + (v.n) [DJ =J, (66 b) 
sena [H,] + (v.n) [D,] = us (66 c) 
cosa [D,] + sena [D,] =p, (66 d) 


As equações (66 a) a (66c) equivalem em nota- 
ção vectorial a: 


nX [HH] + (v.n) [D] =J, (67 a) 
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A equação (66d) a 


n. [DI =p, (67 Db) 


Para as condições fronteiras 
n'X [E] =0 
ns. [B] =0 


obtinham-se por processo idêntico as condições fron- 
teiras (56 b) e (56d). 
Somos assim conduzidos às condições frontei- 
ras (56) mas agora procedendo com generalidade. 
Conclui-se assim que apenas intervém a compo- 
nente normal da velocidade v. 


(v.n) 


Deste modo, para um referencial em que as 
descontinuidades electrodinâmicas se confinem a uma 
superfície constante no tempo: 


F(r)=0 


as condições fronteiras para outro referencial que 
se desloque tangencialmente sobre essa superficie 
são formalmente idênticas âquelas que se observam 
para os corpos em repouso. 

Deste modo também fica clarificado que apenas 
tem interesse prático a consideração da velocidade 
normal 


9F' 
Va a (v.n)n 
IV F| 


E também fica esclarecido que a superfície 
F(r,t)=0 


nada tem a ver com o deslocamento da matéria, 
haverá casos em que coincide, noutros não. 

F(r,t) = 0 é uma superficie que descreve, no 
tempo do laboratório, o lugar geométrico das descon- 
tinuidades electrodinâmicas. 

É assim uma superfície abstracta. 

Este ponto que é importante é salientado por 
Costen e Adamson no trabalho referido na Biblio- 
grafia. 


t. CONCLUSÕES 


— Demonstrou-se que as condições fronteiras 
numa superfície em movimento se podem obter a 
partir das condições fronteiras dos corpos em repouso 
utilizando a Transformação de Lorentz que imobi- 
liza num determinado referencial S' (r',t') o ponto 
P (rt), em movimento, da referida superficie. 

— As condições fronteiras assim estabelecidas 
coincidem com as condições fronteiras obtidas a 


partir da forma integral das equações de Maxwell, 
utilizando os teoremas de Helmholtz associados a 
superfícies e contornos deformáveis no tempo e seus 
corolários. 

— Nas referidas condições fronteiras apenas 
intervém a velocidade normal à superfície 


-— Para o estabelecimento das referidas condi- 
ções fronteiras houve que estabelecer as leis de 
transformação das densidades de carga e de corrente 
superficiais. 

— Com generalidade, estabelecem-se as seguin- 
tes leis de transformação: 


; | cosa” 
J sv OT DM E cos q Esé 
PR sen q 
Ra a A Pa (68 a) 
; J cosa" = sena 
FT cosa — Jstrsma  (68b) 
ca): cos q 
E =y laçe —|Jy—— = 
: Les c? cos q 
V.d, sen q 
=y [p— ——— (68 c) 
Les c? sen q 
W=<L(in,v) u=<L(n,v) (68d) 
nd. =0 (68 e) 


Y Cosa'/cosa=sena';/seng-— tana'= ytana (68 T) 


As condições fronteiras num referencial que se 
mova tangencialmente em relação a uma superfície 
fixa são formalmente idênticas às dos corpos em 
repouso, 

— A superficie F' (r,t) sede de descontinuidades 
é uma superfície abstracta lugar geométrico, no 
«espaço», das descontinuidades que se observam no 
tempo «t» do laboratório. 


— Finalmente convém salientar que se postulou 
que os integrais 


[as v dp! 
| me a —— —— x'dx' 
" or e 5) A 


O 


tendem para zero quando N > O 

Fisicamente o facto é aceitável mas convirá, 
caso a caso, investigar o modelo que permita micros- 
copicamente construir as densidades de corrente e 
de carga e bem assim as leis de variação expectáveis 
no espaço e no tempo de J' e p' 

É matéria cuja análise não incluimos no pre- 
sente trabalho. 

— Não queremos deixar de salientar que, no 
sentido de Sommerfeld, a velocidade 


(v.njmn 


pode ser variável no espaço e no tempo. 

— Finalmente, chama-se a atenção para o facto 
de as condições fronteiras serem um dado importante 
para resolver com elegância determinados problemas. 
Procedendo assim hã que juntar às condições fron- 
teiras as relações constitutivas que se observam no 
referencial do laboratório. 
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Alguns aspectos da correlação de 3.º ordem 
e do bispectro dos processos estacionários 


F. E. REBELO SIMÕES * 


RESUMO 


Introduzem-se, como generalização da 2.º ordem, os conceitos de correlação de 3.º ordem e de espectro 
de 3.º ordem (bispectro) dos processos estacionários, estudando-se as suas simetrias, Faz-se em seguida 
idêntica extensão para correlações e espectros cruzados. Estabelecem-se as relações entrada-saida de siste- 
mas lineares invariantes para correlações e espectros. Considera-se o caso particular dos processos gaussianos 
estacionários de média nula, para os quais se verifica que o bispeciro é nulo, deduzindo-se dai as aplicações 
práticas desta análise. Interpreta-se a seguir o bispectro como resultado de filtragens lineares, à semelhança 
do que se faz para o espectro de potência. Finalmente indicam-se alguns tópicos sobre a computação do 
espectro de potência e do bispectro, comentando-se os prós e os contras da análise espectral de 3.º ordem, 


SUMMARY 


As generalization of 2"! order analysis, the concepts of 3'“ order correlation and 3'* order spectrum 
(bispectrum) of stationary processes are introduced, and their symmetry properties are studied. Identical 
extension for cross correlations and cross spectra is made. The input-output relations for correlations and 
spectra are established in linear time-invariant systems. The gaussian case is considered in which the bis- 
pectrum is shown to be zero. Some practical applications of this fact are presented. Next the bispectrum 
is interpreted as a result of linear filtering operations, like power spectrum. Finalh; some topics on the 
computation of power spectrum and bispectrum are indicated, with comments on the advantages and limita- 
tions of 3'“* order spectral analysis. 


0. INTRODUÇÃO 


Os processos estocásticos aparecem na Engenharia como modelos matemáticos adequados de certas 
grandezas aleatórias dependentes de um ou mais parâmetros. No presente artigo consideraremos apenas um 
parâmetro contínuo (t) que admitiremos ser o tempo e restringir-nos-emos a grandezas reais. 

Assim, uma amostra ou realização do processo terá o aspecto de uma função aleatória real x(t), que 
é costume designar por sinal, principalmente na linguagem das Comunicações, 

Notemos desde já o carácter duplamente aleatório deste sinal x(t), no tempo e no conjunto. Por um 
lado x(t) varia em cada instante t de amostra para amostra de forma aleatória, constituindo uma popu- 
lação x (t), x,(t),... x (t),... com determinadas propriedades estatísticas que deverão ser conhecidas para 
que o processo possa ser analisado. Por outro lado x(t) varia em cada amostra x(t) de instante para 
instante de forma — em geral(1) —- igualmente aleatória, obedecendo os valores x(t,), x(t,)... x(t )..., a 
certas leis que precisarão também de ser conhecidas. 

Este duplo comportamento tem de ser visto de forma global, pois os dois aspectos não são inde- 
pendentes. 


* Centro de Análise e Processamento de Sinais das Universidades de Lisboa, e Instituto Superior Técnico. 
() Este caso geral aqui considerado corresponde ao chamado processo puramente aleatório. 
Manuscrito recebido para publicação em 11/4/77. 
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Assim, por exemplo, o conceito de estacionaridade, de 2.º ou de 3.º ordens, diz respeito ao modo como 
evolui no tempo a média de conjunto respectivamente dos produtos de 2.º ordem x(t) x(t + +), com « fixo, 
ou de 3.º ordem x(t) x(t+<,)) x(t+«<,), com Tt; € T, fixos. A estacionaridade de qualquer das 
ordens existirá quando a respectiva média acima referida não depender de t ('). 

Por sua vez o conceito de ergodicidade (para o caso estacionário) tem a ver com a comparação 
entre os valores médios de conjunto atrás definidos e os valores médios tomados ao longo do tempo para 
as diferentes realizações. A existência de ergodicidade implicará que estas médias temporais devem ser 
iguais para todas as amostras e coincidir com as médias de conjunto correspondentes. Por exemplo, se um 
processo estacionário de 2.º ordem é ergódico tem-se necessariamente 


T/2 
(2) 1 E 
E x(t) x(t++) = lim— |xít) x(t+4 +) dt (0.1) 
T-—+>0c0 T 
—T/2 


No que vai seguir-se consideraremos apenas processos estacionários até à 3.º ordem, (em sentido lato) 
isto é, processos para os quais se tem, para qualquer valor de t,, 


Ex()=Ex(t+t)=m, (0.2) 
E x(t) x(t +) =E x(t+t) x(ttHtitr)—m, lr) (0.3) 
Ex(x(t+tz)x(ttr,) =Ex(t+t) x(t+tt+z) xt+Httet)=m, (moto) (0.4) 


Estas médias não dependem, portanto, da origem do tempo, sendo constantes para a 1.º ordem e 
dependentes apenas dos intervalos entre os instantes considerados, para as 2.º e 3.º ordens. 

Nas expressões apresentadas m, é o momento de 1.º ordem ou valor médio, m,(t) é a função momento 
de 2.º ordem e m, (t,,7,) à função momento de 3.º ordem. Se nas expressões (0.3) e (0.4) subtrairmos a 
cada factor o valor médio m,, obteremos as chamadas funções momentos centrais da ordem respectiva. 

Em particular suporemos no que segue m =, o que torna as funções momentos coincidentes com 
as funções momentos centrais. Estas funções coincidem ainda com as chamadas funções de correlação (ou 
simplesmente correlações). Acima da 3.º ordem o conceito de correlação difere do de momento central, mas 
esse caso não será aqui abordado (?). 

No presente artigo consideraremos ainda ($ 8.) o caso dos processos ergódicos até à 3.º ordem. Isto 
implicará, no caso m, = 0, as seguintes igualdades válidas para todas as realizações do processo [cf. (0.1)]. 


T/2 
EXMO)= lim É ) x(t) dt=— 0 (0.5) 
T-+>00 T 
—T'/2 
T/2 
Read 
E x(t) x(t + +) = lim = [x(t) x(t +) dt= q (T) (0.6) 
Roo Do 
—T/2 
T/2 
E x(t) x(t+—,) x(t+ +,) =. lim a) x(t) x(t ++) x(t+t,) dt=y (Tot) (0.7) 
T>o T E 
—T/2 


em que designamos por q (7) e y(z,,T,), para evitar os índices 2 e 3, as correlações de 2.º e 3.º ordens. 

As igualdades anteriores significam que o valor médio (neste caso, nulo) e as correlações de 2.º e 
3." ordens se podem obter (teoricamente) através de operações — médias temporais — efectuadas a partir 
duma única amostra, suposta de duração infinita. 


(1) Este tipo de estacionaridade é por vezes chamado fraco ou em sentido lato. 

(2 E fit) g(t) significa «valor médio de conjunto do produto f(t) g(t)». Evitaremos usar parêntesis a seguir ao gim- 
bolo E para simplificar a notação, 

(1) Ver por exemplo Stratonovich, 1963, Cap. 1. 
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A hipótese de estacionaridade é fundamental para a definição dos conceitos de espectro de frequência 
de 2.º e 3.º ordens e para o estabelecimento das propriedades aqui tratadas das correlações e respectivos 
espectros. 

Quanto à hipótese de ergodicidade, ela permiti--nos-á abordar o problema da obtenção experimental 
dos espectros de 2." e 3.º ordens por meio de operações realizadas sobre amostras isoladas. 


1. CORRELAÇÃO DE 2.º ORDEM . SIMETRIAS 
Como vimos, tem-se 
q (T) = E x(t) x(t + q), com Ex(t) = 0 
Usando em (0.3) t, = —r, resulta 


(it) =E x(t—r) x(t) = q (—<) (1.1) 


isto é, a correlação de 2.º ordem é uma função par. (Notemos que estamos apenas a tratar de processos 
reais. Caso contrário, o conceito de correlação obrigaria a uma ligeira alteração da definição de correlação, 
que passaria a ser, por exemplo, q (7) = Ex*(t) x(t ++), em que x*(t) representa o conjugado de x(t). 
Mas não consideraremos aqui o caso complexo). 


2. ESPECTRO DE 2.º ORDEM OU ESPECTRO DE POTÊNCIA 


Na maior parte dos processos puramente estocásticos (sem componentes periódicas) que servem de 
modelo aos sinais encontrados na Física e na Engenharia, verifica-se que 


lim q(r)=0 
5:65 (2.1) 
de tal modo que se tem 
Do 
[ |p(r) [dr < 00 (2.2) 


Esta condição é suficiente, como se sabe, para a existência de transformada de Fourler de q (T) 


& (w) = fo cas ga (2.3) 
tendo-se 


o (x) = — + (mw) e dy (2.4) 
ke 


(1) Tratamos aqui da 2.º ordem apenas como meio de melhor compreender, por extensão natural, as propriedades das 
grandezas de 3.º ordem que são objecto deste trabalho. 
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A grandeza + («y) é, por definição, o espectro de 2.º ordem do processo, e tem o significado duma 
densidade espectral do valor quadrático médio E xz(t), usualmente chamado potência ('), como resulta de 
se ter 


e 


E x:(t) = q (0) = —— | (mo) du (2.5) 


. 


É evidente que da propriedade (1.1) de simetria de w (+) resulta idêntica propriedade para  (ím) 
bi—mw)=*(m) (2.6) 


O espectro de potência é uma função real e par. Pode provar-se ainda que é não-negativa (7), como 
seria aliás de esperar do seu significado físico de densidade espectral de potência. 


3. CORRELAÇÃO DE 3.º ORDEM. SIMETRIAS 


A definição de correlação de 3.º ordem e a hipótese de estacionaridade vão permitir-nos obter as 
propriedades de simetria desta função. 


Com efeito, de (0.4) resulta, por conveniente escolha de t , o conjunto de identidades 


pls ta) = Êta T) 

= Gasto Ta) 
fra tata) (3.1) 
A e 


o Al Cana OS Ra O 
como se pode verificar, fazendo por exemplo t = —+ em (0.4). 
Tem-se então 
pit tT,) = Extt-z)x(t)x(t-t +t,)=y(—t,tT.—T,), que é a 2." das identidades (3.1) 


Estas relações (3.1) exprimem afinal a equivalência entre as expressões que se obtêm para a correlação 

de 3." ordem y (tT,, 7.) escrita sob a forma 
pit to) =E x(L)x(t) x(t) =p (t—t o, t=t,) (3.2) 
quando se permutam os instantes t,, k,, t. das 3! maneiras possíveis, alterando apenas a ordem dos 3 factores. 
t, virá 


Por exemplo, se tomarmos a permutação CU, U ta=t,t ! 


y (To To) = Ext, x(t)) x(t) — Ex(t) x(t,) x (ti) — 


tt tt) =y(t—t ,t-t)=y(z,-t,,—+,), que é a 3.º das identidades (3.1) 


Esta forma de proceder é aliás a mais aconselhável para ordens superiores à 3.º 


As simetrias de y (+, ,7t,) dadas por (3.1) correspondem à partição do plano (t, T,) em seis sectores 
equivalentes, definidos pelos eixos principais e pelo bissector dos quadrantes ímpares. (Fig. 3.1). 


(1) Se x(t) é uma tensão aos terminais duma resistência, xt) é proporcional à potência instantânea. Nos processos 
ergódicos E xt) será proporcional à potência (média) do sinal. Por isso é usual chamar a P(w) o espectro de potência. 
0º) Ver por exemplo Davenport & Root, 1958, Cap. 6. 
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L, 


A pi T.) he, 


P, (TT, Ly) 


Fig. 3.1 


Com efeito, se o ponto E, (x, + Ta) percorrer o sector 1 (1,2 1,2 0), por exemplo, o ponto homólogo 


e 


Pitta E tu —tT,) que resulta da 3.º identidade (3.1) percorrerá o quadrante 4, por se ter a 
implicação 


>,>20=> 7,20,7,4<0 
Da mesma forma se obtêm os restantes sectores. A cada ponto dum sector corresponde um ponto 
homólogo em cada um dos outros, no qual a correlação | (r,, To) toma o mesmo valor, Deste modo, o 
conhecimento de Ee; Rg) num dos sectores é suficiente para determinar esta função em todo o plano. 


4. ESPECTRO DE 3.º ORDEM OU BISPECTRO 


Tal como para a 2.º ordem verifica-se, em grande parte dos processos encontrados na prática, a 


condição 
(e 
ff Wit, sto) | dr, dr, < 00 (4,1) 
a 


que é suficiente para a existência de uma transformada dupla de Fourier de y (To To): Tem-se assim 


Do 


já ES td TT To) 
y (19, 4) = | BAT Ta) 8 dr, dr, (4.2) 


1 (mt, + w, To) 
= | | Vlw, t,) e do, de, (4.3) 
— 00 


bit Ta) 
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A grandeza W (mw, ,(,) é, por definição, o espectro de 3.º ordem ou bispectro do processo e parece 
ter o significado duma densidade espectral a 2 dimensões do valor cúbico médio E xs(t), como resulta de 
se ter 


So 

1 

Ex'(t)=y (0,0) = =) | E (mw, to,) do, de, (4.4) 
— O 


Na realidade o bispectro W («),, ty ,), dO contrário do espectro de potência & (ww), é em geral complexo, 
e apenas a sua parte real contribui efectivamente para o valor médio. Veremos isto adiante. 

Para uma interpretação da parte imaginária, ver o 58. 

As simetrias do bispectro W («,, tu 53 pdem obter-se como resultado das simetrias da correlação de 


3.* ordem y(T,.T,): 


Assim, por exemplo, da relação de simetria 


Pratt) =YtT-t, + —t4) 
tem-se 


co oo 
=] (td; T; + toTo) | =] Ltd, Ty *F Wa Ta) 
E (ww, tw,) = [fo (T,,7;) O dr, dr, -| PCT Ts Tt) O dr, dr, 
“ow n 
— (E) 


mm (JC) 
Fazendo t-t=T, e —t,= 7', O último integral (1) pode escrever-se 
co 
À 8 Eros to Cio, + toh Tr] 
do PARTO) O dr", dr, = F (m,, —t, —to,) 


De forma semelhante se podem obter as restantes relações, resultando o conjunto de identidades 


E (to,.0,) = E (mm, st) 
=Y(wm, —tw, —t,) 
=+(—w —Wootd) (4.5) 
== (tdo 04 0) 
=F(—w —tW,o t) 
Pode provar-se (2) uma regra simples para escrever imediatamente estas equivalências: consiste em 
interpretar W (tm, . w,) como possuindo 3 variáveis «w, to, . td, ligadas pela relação w, + w,+tw,=0, e pro- 


ceder à permutação destas 3 variáveis, explicitando como argumento de Y apenas as 2 primeiras (apenas 
2 são independentes), e substituindo sempre w, por =, —Ua 


Por exemplo 


Elo) <> ww) = Tlm, mat) <LD> 
<> tu) = (mw —to» to) 


que é a 5.* das identidades (4.5). 


Fi) (+, ' r)) 
() Notar que o jacoblano ————— & igual a 1. 
Er, à Tr") 


(2) Ver por exemplo Brillinger & Rosenblatt, 1967. 
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Estas relações de simetria (4.5) correspondem à definição de seis sectores numerados de 1 a 6 
no plano (wo, w,) (Fig. 4.1). 


W, =, 


SÉ 7 
SN 


Fig. 4.1 


Com efeito suponhamos um ponto P (w,,«w,) que percorre o sector 1 (w, 2 tw, 2 0). A 5.º identidade 
(4.5) por exemplo, define um ponto homólogo P. (w',, w',) = P.(—-w,—tw,, t),), onde o bispectro toma o mes- 
mo valor que em P, e que percorrerá o sector 3, dada a implicação 


1 
W, 2 w,2 0 =» 0<w,<— wa 


Da mesma forma se obtêm os sectores 2, 4, 5, 6.. Estes seis sectores não preenchem, no entanto, 
completamente o plano («w, w,), restando outros tantos verticalmente opostos aos anteriores, onde o bis- 
pectro se pode calcular facilmente (sectores indicados com asterisco). Isto deve-se à existência da relação 
de simetria óbvia 


E (wo tw) = W* (6, q (0,) (4.6) 
que resulta directamente da definição (4.2) de W (mw, , to,). 


Assim 
co 


JC, T+ mw, To) 
F(—w, —w,) = EN) Y lt to) & dr, dr, = ºC,» o) 


— (o) 


dado o carácter real da grandeza y (T,,T,). 


Nestas condições, o conhecimento do bispectro num destes doze sectores determinará completamente 
o seu conhecimento em todo o plano («w, w,). 

A Fig. 4.2 mostra como se podem obter os cinco pontos homólogos de cada ponto P do sector 1. 
Os seis pontos restantes (onde o bispectro toma um valor conjugado do anterior) obtêm-se por simetria 
em relação à origem. 
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no tm 


ty, 


Fig. 4.º 


Em face do que fica dito podemos escrever a relação (4.4) utilizando apenas metade do plano (ty, w,), 
desde que substituamos Y por 2 Re (W). Por outro lado as simetrias estudadas reduzem-nos o domínio de 
integração a um dos seis sectores numerados de 1 a 6, dada a equivalência entre os pontos homólogos dos 
sectores e o facto de o Jacobiano das várias transformações ,,,, tw, — wo tw", entre os sectores ser igual 
à unidade. 


Nestas condições, escolhendo por exemplo o sector 1, vem 


a CO tb, 
Me | | Re (% (9, » t3,)) deo, dto, (4.7) 
wo, =0 o 


que mostra que é a parte real do bispectro que se deve interpretar como a densidade espectral do valor 
cúbico médio, ou seja, a contribuição de cada par de frequências para esse valor. (Ver ainda o $ 8). 


5. CORRELAÇÕES E ESPECTROS CRUZADOS 


Em muitas aplicações, nomeadamente no estudo dos Sistemas, tem interesse a extensão dos conceitos 
anteriores de funções de correlação e respectivos espectros a mais do que um processo, 


Assim, se forem (x(t)) e (y(t) ) dois processos estacionários de 2." ordem, de valores médios nulos, 
define-se como correlação cruzada de 2.º ordem a grandeza 


Palt, E) =E x(t)y (t) (5.1) 
em que x(t) e vít) são realizações dos processos referidos. 


Se esta função Pa lt, st,) apenas depender da diferença t, —t,, diz-se que os dois processos são 
estacionariamente correlacionados (de 2.º ordem). 
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Podemos então escrever aquela função como de uma só variável, tal como fizemos para o caso do 
processo único. 


ottt)=a (D)=ExbDy(t+r) (5.2) 


Nestas ecndições define-se o espectro cruzado Do, (1) 


00 : 
=" | ()T 
Py(w) Lost) e dr (5.3) 
— OO 
tendo-se então 
[9.2] 
1 jwr 
Pay tt) = ato e de (5.4) 


para o que basta a condição, encontrada na prática, 


o 


[ | Pay lt) | de < 00 (5.5) 


= 


— 0 
Quanto às simetrias, repetindo os raciocínios usados nos $ $1.e 2., obtém-se, para a correlação cruzada 
Pyt—T) = Pb) (5.6) 
e para o espectro cruzado 
Sw) = 8, (09) (5.7) 


Da mesma forma, se tivermos três processos ( x(t) ), (yít)) e ( z(t) ) estacionários de 3.º ordem, de valores 
médios nulos, define-se como correlação cruzada de 3.º ordem a grandeza 


Vryatts to t)=E x(t)y(t,) zht,) (5.8) 

Se esta função depender apenas das diferenças entre os instantes considerados, diz-se que os três processos 
são estacionariamente correlacionados (de 3.º ordem). 

Podemos então escrever, tal como no caso de um único processo, e referindo as diferenças de tempo 


ao instante t, que representaremos simplesmente por t, 


Peatt+tro tt) =Ex(Dy(t+r)zt+T,) =yryalTio To) (5.9) 


XYE 


Nestas condições, e verificando-se habitualmente a condição 


oo 
f si | Waya(T1 + Tó) | de, de, < 00 (5.10) 
— DO 


define-se o bispectro cruzado W, (mw, w,) | 


o 
| 5 C08 Td QT) 
Eça MD, + D,) = [ [ta Ta) € dr, dr, (5.11) 


e ÃO) 
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tendo-se então 


LE) 


1 A d Ceu, Tt; + tds Ta) | 
Vaya l Ta p Ta) = Tr | Po tt, e (do) e de, de, (5.12) 


— DO 


Quanto às simetrias, se repetirmos o raciocínio do 8 3., obteremos para as correlações cruzadas o seguinte 
conjunto de identidades 


PrayÃ ta o Ty) 

Fi To TE Ti) 

E faltas) (5.13) 
E fat To tg) 

a Vaxyh— Ta , Ti To) 


Vrya (Tt; ' Ta) 


Procedendo agora para o bispectro cruzado como fizemos no $ 4. obtêm-se as seguintes simetrias 


Fra ttidy ' td) e E cry tdo , to, ) 
E att), » —t) —tda) 
E Eogt—U — tg 004) (5.14) 
a E sxtta » lo — da) 
— EM — Was | (05) 


Resta dizer que todos estes conceitos e propriedades ,tanto na 2.º ordem como na 3.º ordem, degene- 
ram nos anteriormente tratados do processo único, quando se faz x(t) = y(t) — 2.º ordem — ou x(t) = 
= y(t) = z(t) — 3.º ordem. 


6. RELAÇÕES ENTRADA-SAIDA DE SISTEMAS LINEARES INVARIANTES, PARA CORRELAÇÕES E ESPECTROS 


Se for h(t) a resposta dum sistema linear invariante no tempo à função delta de Dirac a(t), sabe- 
mos que se tem 


y(t) = | x(t— ph (u) dy -fa (ujhit—u) du EL) 
= 06 — (O) 


em que x(t) e y(t) são respectivamente os sinais de entrada e de saída do sistema. 

Dado que, como dissemos atrás, as correlações cruzadas e respectivos espectros contêm como caso 
particular as correlações e espectros simples, trataremos aqui do caso mais geral de dois ou três sistemas, 
conforme a ordem a estudar. 

Assim, consideremos dois sistemas lineares invariantes, definidos pelas respostas h (t) e h,(t), e cal- 
culemos as correlações cruzadas de 2.* ordem à saída, supondo os processos de entrada (x (t)) e (x,(t)) 
estacionariamente correlacionados e de médias nulas (Fig. 6.1). 


Fig. 6.1 
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Teremos em primeiro lugar, em virtude de (6.1) 


E y(b=H/(0) Ex (t)=0 


(6.2) 
E y(t=H(0) E x,(t)=0 
em que H (w) e H,(w) são as funções de transferência dos sistemas, dadas por 
[2 a) 
— Ttut 
H(«) = [14 e dt 1=1,2 (6.3) 
— «DE 


As relações (6.2) mostram-nos que os valores médios dos sinais de saída dos sistemas são nulos como os de 
entrada, pelo que as correlações cruzadas à saída serão dadas por 


o 
Oy trt) = Eyilty) Volto) = ef [x (t=) x, (t>v) bh, (pu), (v) dy dy 
1W2 
cms 
o 
(ão 6.4 
== E (t-uot,—v) h (gu) h,(v) dy dy (6.4) 
E 
00 


Vemos que 2, yo t,) dependerá de t, e t, através da diferença t,—t,, se e. Ro —t o, t,—v) também 


o fizer. Isto é, a cotndicsiidiado correlacionada dos processos de entrada Pa a dos processos de 
saída (2). 

No caso presente suporemos (x,(t)) e (x,(t)) estacionariamente correlacionados. Nestas condições, 
designando Raça por 7 e v—| por É, a integração (6.4) fornece a relação 


co 
= a À (£) (6.5) 
Ca (+) for (=) Ph, E) dg 
— 0 
em que 
[88] 
Ph, (D=|h()h,(u+v) dy (6.6) 
Es 
Conclusão: 


A correlação cruzada de 2.º ordem dos processos de saída obtém-se da correlação correspondente de 
entrada, através da sua convolução com a função de correlação cruzada das respostas h (t) e h,(t) dos 
dois sistemas. 


(') Prova-se que se se verificar a condição p E |x(t) | dt «co para um processo ( x(t) ), se pode escrever 
a 


E P x(t) dt = fE x(t) dt, tendo-se ainda f| =(t) | dt, < 60. Ver Davenport & Root, 1958, cap. 4, pág. 65 
a 


(2?) Esta condição é suficiente, mas não necessária, 
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Estes produtos de convolução têm uma correspondência no domínio da frequência em termos de pro- 


dutos simples, ou seja 


+ (u=P (oe (tm) 
Fy7a EE, hyh, 


em que ETA (tw) é a transformada de Fourier de Ph (£), dada por 


Ta toa 


(19) =H (—w) H,(m) 
4h, 


= H *(mw) H,(w), por ser h(t) real 


Conclui-se assim que o espectro cruzado à saida dos dois sistemas é igual ao 


(6.7) 


(6.8) 


produto do espectro 


cruzado das entradas pela função H (—.) H,(ww), que podemos denominar por função de transferência cru- 


zada de 2.º ordem. 


Notemos que se podem particularizar estes resultados para os casos h =h, e/ou X =X, 
A generalização para a 3.º ordem destas relações entrada-saíida não otérece dificuldade. 
Assim, considerando três sistemas lineares invariantes (Fig. 6.2) definidos por h (t), h,(t) e h,(t) e 


três processos de entrada (x (t) ), (x,(t)), (x, (t)) 


Fig. 6.º 


estacionariamente correlacionados de 3.º ordem, e de médias nulas, tem-se 


4 (Lot t)=E y(tO)yÃt) vilt,) 
FF 253 


(6.9) 


e por um raciocínio semelhante ao anterior, concluiremos primeiro da estacionaridade correlacionada dos 


processos de saída, obtendo finalmente a relação 


[6 
ovas a Ta) E | | E CEM Pan A 8a) de, de, 
— 00 
em que 
Do 
Phu Ea ) = ) h, (pt) h, (q + E)h, (u + E.) du 
— co 
ou seja: 
26 


(6.10) 


(6.11) 
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A correlação cruzada de 3.º ordem dos processos de saida obtém-se da correlação correspondente 
de entrada através da dupla convolução desta com a função de correlação cruzada das respostas h(t), 
h,(t) e h.(t) dos três sistemas. 


Tal como para a 2.º ordem, estas relações traduzem-se no domínio da frequência por simples produtos 


(ro, ' Wa) 


qr (ww, ,W,) =Y (o, to) F 
Uia o Vig 2 hhh, (6.12) 


P4F95's K,X,X, 


em que w giinis (o, » w,) é a transformada dupla de Fourier de E mn (£, + $,) dada por 
| 


IOLTa 1 
Y (ow) =H(—w,—w,) H, (w,) EH (w,) (6.13) 
hjhh, 1 2 1 1 ia é 1 3 2 
=H*(w, +w,) How) H(w,), por ser h (t) real 


Também estes resultados são susceptíveis de particularização, tanto no que se refere à igualdade entre 
os sistemas (todos ou apenas dois deles), como no que se refere aos processos de entrada (x =x, =X, 
ou, por exemplo, =X, 75X%,) 


7. O CASO GAUSSIANO. IMPORTÂNCIA DA 3.º ORDEM 


Como se sabe, um processo estocástico diz-se gaussiano se, para qualquer conjunto finito de instan- 
tes (Loto t, ), as variáveis aleatórias (x,,X,,...X,), com x =x(t ), têm uma distribuição con- 
junta gaussiana ou multinormail, 

Isto significa que se tem para a função densidade de probabilidade conjunta a expressão (1) 


N N 
1 | 
P(X,X, ex) =—"——— > exp)  — > b3 Yu —mo (—my) (7.1) 
N/2 1/2 o a E 
(25) (detA) p=] k=l 
em que 

m, = Ex, (7.2) 

A é a matriz das correlações, de elementos À, dados por 
A = E(x—-mp) (x—m) (7.3) 


e vY,, São os elementos da matriz inversa de A, [= AT, 

Embora neste trabalho nos restrinjamos aos processos estacionários de média nula, a análise que se 
segue é válida no caso geral, 

Vamos mostrar que a correlação de 3.º ordem é nula nestes processos. 

Para isso (ver Apêndice), basta fazer uma mudanca de variáveis, através da transformação linear 


ho 
= mem ee Dao N (7.4) 
k==1 


em que fizemos 


JH = %4—m, (7.5) 


Esta transformação, por ser linear, mantém o carácter multinormal das novas variáveis. (Ver A 19). 
Se designarmos por A a matriz, de elementos a da transformação (7.4), a nova matriz das corre- 


lações das variáveis y', é dada por 


lk * 


T 
A =ALA (7.6, AT, 


(1) O det A é sempre positivo (ver (A 28), Apêndice). 
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É possível diagonalizar esta matriz através duma escolha apropriada de A, obtendo-se para as novas 
variáveis y', a distribuição 


N 
PY IV y= HH P(',) (7.7, A 25) 
j=1 
em que 
1 É 
PO) CEA Pp (—y',/2h,) 


sendo à, os elementos da matriz diagonal A', dados por à, = Ey*,º 


As novas variáveis y', são portanto independentes, o que nos vai permitir calcular facilmente E y y,y, 
e verificar que a função de correlação de 3." ordem dum processo gaussiano, mesmo não estacionário, é 
identicamente nula. 

Com efeito, considerando 3 instantes quaisquer t,, t,, t, e designando, como já o fizemos, as variáveis 


= 


x(t)— Ex(t),x(t,)— Ex(t,) e x(t)—-Ex(t,) POr 9,)344%,+ É possível obter por uma transformação 
conveniente, de acordo com o que ficou exposto, 3 novas variáveis y',,Y',,Y',,» obedecendo a (7.7, 7.8), 
com N — 3. 


A correlação de 3." ordem procurada será dada por 


A om 


v (ts tos to) = By, Ja Va E | | | LÊ Ja Vs P(y, For SD) dy, dy, dy, (7.9) 


—— DO 


Dado que os y, também se exprimem nos y', através duma transformação linear (expressa pela ma- 


-1 T ' 
triz inversa de A, A = A, ver Apêndice) o produto y,y.,Y, é constituído por parcelas que, à parte coefi- 
cientes constantes, são de forma y',3, y,2 Y, (1Ak) e y', Y', Y',- Todas estas parcelas vão ter contribuição 
nula para o cálculo de BY Fe 


Com efeito, em virtude de (7.7) tem-se, por exemplo 


oo 
E y',*= ff y'3 Pl) PCS) P(y',) dy”, dy”, dy', 
—o 


o So o 


= t y 3 P(Y,) dy, ) p(y,) dy, f p(y,) dy”, = 0 (710) 


LO 


— 00 — 00 —00 


por ser nulo o 1.º integral (integranda ímpar). 


Da mesma forma 


oo oo co 
E y,2y,= f YºPOyp) dy, | Ply',) dy, À Fay) d7,=0 (711) 
por ser agora o 3.º integral que se anula. 
Finalmente 
oo 
3 
Ry,yg= T EZEIcA dy, =0 (7.12) 
Ll= 1 a 
— OO) 


por serem os três factores nulos. 
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Está assim provado, como pretendíamos, que a correlação ou momento central de 3.º ordem dum 
processo gaussiano, não. necessariamente estacionário, é nula. 

Este mesmo raciocínio aplicado a ordens superiores à 3.º mostra que todos os momentos centrais de 
ordem impar são igualmente nulos. Acima da 3.º ordem as correlações diferem dos moments centrais, sendo 
definidas doutro modo. Pode provar-se que todas as correlações a partir da 3.º ordem são nulas para o 
processo gaussiano. (Este anulamento das correlações é aliás uma forma axiomática de definir o processo 
gaussiano, utilizando a função característica) (1). 

Só para a caso estacionário se pode fazer a descrição dos processos em termos da frequência. É neste 
contexto que se define o bispectro, como se viu (5 4.). Podemos pois concluir que o bispectro dos processos 
gaussianos estacionários é nulo. 

O interesse da análise de 3.º ordem, nomeadamente da análise espectral, hoje mais exequível por via 
computacional do que a análise de correlação (dada a existência de algoritmos muito rápidos de análise 
de Fourier como a FFT — Fast Fourier Transform) reside precisamente na sensibilidade daquelas gran- 
dezas de 3.º ordem relativamente ao processo gaussiano, para o qual se anulam. 

Isto permite avançar em 2 tipos de problemas: 


1) Detectar a presença dum processo não gaussiano num fenómeno qualquer em estudo (2). 
2) Detectar a natureza não linear dum fenómeno, admitindo que as suas causas possam satisfazer 
uma lei normal. 


Este 2.º tipo de problemas é o que aparece com mais frequência, principalmente em Geofísica (3). 


3. O ESPECTRO DE POTÊNCIA E O BISPECTRO COMO RESULTADOS DE FILTRAGENS 


Vejamos como se pode interpretar o espectro de potência e o bispectro dum processo estacionário em 
termos de filtragens lineares de banda estreita. 

Suponhamos então que fazemos passar o sinal x(t) por dois sistemas lineares invariantes constituindo 
filtros passa-banda muito estreitos, centrados nas frequências ww, e w, Vamos admitir que estes filtros são 
ideais, de largura A, definidos por 


H(w) — ela w elw,—4/2, w, +4/2] (8.1) 
veto,  wel-w—A/2, —w,+4/2] 
= 0 fora destes intervalos 


com 4, constante 


H.(w) supõe-se análgo a H (mw), substituindo w, porw, e 9, por 8,. 
Estes filtros, embora reais (8.1), não são evidentemente realizáveis fisicamente, apenas se podendo 
realizar matematicamente. (Fig. 8.1). 


Fig. 8.1 


(1) Ver Stratonovich, 1963, Cap. 3. 
(:) Ver Rebelo Simões, 1971. 
(1) Ver Hasselman et al, 1963; Mac Donald, 1965; Madden, 1964; Haubrich, 1965. 
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Calculemos agora o valor médio do produto y (t)y,(t) das saídas destes filtros, admitindo que 
os intervalos definidos em (8.1) correspondentes aos dois filtros não são respectivamente disjuntos (Fig. 8.2) 


Ko. 
-Wg (8) | 
MA. 
+= 4 = e | 
do O had 5 


Podemos exprimir E y,(t) v.(t) em função do espectro de potência cruzado & ae (w), e este, por 
já 12 
sua vez, no espectro de potência de entrada & (mw). 


Com efeito, atendendo às relações (5.4), (6.7), (6.8), tem-se 


Do 
(1) 1 
RM Fa a | ga a Coe 
— co 
So 
= | tw) H(-w) H,lw) do (8.2) 
e 
==) 


Pelo teorema da média, poderemos escrever, supondo Wy 2 Wa 


to, + A/Z 
» 
E y,(t) y,(t) — ao db (w,). 2Re | / H(—w) H(w) de (8.3) 
w,—4/2 
em que ww é um valor do intervalo da integranda. 
Calculando este integral, obtém-se 
| wW,— to, + à 
DIA ÃO —=————. cos (a, —9,) & (mw) (8.4) 
q 


Suponhamos que o espectro em estudo P (mw) varia muito pouco no intervalo 4. 


Se fizermos agora varia ty, & 6, admitindo ww, € 6, fixos, a expressão anterior é máxima para ty = 


da ' 
0, = 8, vindo 


â 
TAIT) — E (mw) 


au 


(8.5) 


Nestas condições os filtros são iguais e y,(t) =y,(t), pelo que a operação é na realidade a determinação 
do valor quadrático médio de um sinal filtrado por um único filtro, centrado em w= tw, de largura À e de 
fase 4, constante na banda, tendo-se como resultado que (8.5) se pode escrever 


â 
1 at) = Dl) 
di x “Mg 


ou seja 


+ (uv) = e E y2(t) 
(8.6) 


() E y(O=E y(1)=0, mesmo que Ex(t)-£0, por se ter H,(0)=H,(0)=0 (cf. 6.2). 
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q 
por 


Um novo &s 


pa é 


ELLTIGA. 


X 


TECNICA 


TECNICA X 


a RSRS RO SE o. «io RD = RR 


EstirAL O” 


revista 
internacional 
de luminotecnia 


Publicação trimestral, a córes, com informação 
sobre os mais importantes dominios de 
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Se suposermos que o processo em estudo é ergódico, o valor médio de conjunto referido pode ser 
substituído por uma média temporal, tomada ao longo dum tempo infinito, ou seja 


Sw) = es y2(t) (8.7) 


Na prática podemos substituir estes métodos puramente matemáticos por métodos físicos, para o que 
basta substituir os filtros ideais apresentados por aproximações realizáveis em que 


H(w) = ed 9 (0) nos intervalos referidos (8.8) 
= O fora desses intervalos | 


tendo-se, como anteriormente H (—w) = H*(ww) para assegurar o carácter real de hit) e a fase stw) 
variando adequadamente com w para assegurar a realizabilidade física, 

Nestas condições, e pondo w, = tw,» Gpltw) = 6, (ww), obtém-se igualmente a expressão (8.7). 

E este o principio de funcionamento dos analisadores analógicos de espectro de potência, 

Deve notar-se que além do erro introduzido pela inexactidão do filtro empregado em relação ao ideal, 
temos o erro mais importante resultante da duração limitada do intervalo de avaliação da média. 

Esta limitação é mais acentuada nos métodos analógicos que operam em cada instante sobre um certo 
passado, do que nos métodos digitais, utilizando sistemas fisicamente irrealizáve;s como q descrito .nicial- 
mente, já que nestes métodos se opera em cada instante sobre todo o registo. Mesmo assim, o erro é ainda 
razoável, podendo ser substancialmente diminuído com o emprego de métodos de estimação estatística, 
que referiremos mais adiante a título de informação. 

Vamos agora ver como se pode estender à 3.º ordem a interpretação do espectro em termos de fil- 
tragens de banda estreita análogas às anteriores, mas utilizando três sistemas. 

Para isso consideremos três filtros ideais defin dos como em (8.1) mas centrados respectivamente nas 
frequências w,, to, tw,» € de fases 9,, 0,10, (Fig. 8.3). 


Assim, e nos intervalos correspondentes, de largura 4, será 


8 

H(w) = 03% 
jo - 
H(w=e * (8.9) 

ê 

H,(w) = Pd 


anulando-se H (tw) fora desses intervalos, e tendo-se H,(—w) = H,*(tmw) 


PAGIARIAS 


Fry. 4.3 
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Calculemos agora o valor médio do produto y (t) y,(t) y. (t) das saídas destes filtros. | 


Podemos exprimir este valor em função do bispectro cruzado W (to, » tu,), é este por seu turno no | 
V$F3Y 
bispectro da entrada W (tÔ,, tw). e 


Com efeito, atendendo às relações (5.12), (6.12), (6.13), tem-se (1) 


co 
1 os dire 
E y (t) att) At) — 4 (0,0) = | jr (ty st) ) dis dty ee 
1 Va Va Es (25)? J | psy ab, 2 
— (NC) 
es 
1 ra | 
= a: | | tutto va) Hit, —00,) Holio,) El) des, du, (8.10) 
so 


Devido às condições (8.9) o produto H (w,) H.(m,w,) vai delimitar o domínio de integração aos 4 qua- 
drados figurados (Fig. 8.4), em que se supós, sem perda de generalidade, ww > w, > O 


| 
Fig. 8.4 


Por outro lado, e admitindo tw, € w, fixos e 1, variável, o factor H(—y, — (4) só não se anulará se as 
faixas 


w,—8/2<w, tw, < tw, + 4/2 (8.11) 
=, 4/20, +u<—w + 4/2 


intercectarem os quadrados centrados em (tw ,w,) e (—w.,—tw,) — caso em que « está na vizinhança 
de (y 4 ty, OU os quadrados centrados em (mw. ,—w,) € (—w., tu) — caso em que ty, está na vizinhança 
de ww, — tw. 


E 


() Também aqui se tem Ey (1)=Ey (t)=Ey (=0, mesmo para Ex(t)£0, por se ter H/(0)-=H.(0)=H,(0)=0 
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Consideremos então a situação representada (Fig. 8.5), que corresponde à w = tw, + tw, 


W, 


B 
| Tuscussucen/2 


a AR cmi à | 


5 + w,=-U,+8/2 is 


é 2 


[ 
| w+w,=-w-A/2 


Fig. 8.5 


Se suposermos que 4 é suficientemente pequeno para que o bispectro varie muito pouco na região 
de integração R (zonas tracejadas), o integral (8.10) pode escrever-se 


E (mt) | | H(—mw, —to,) Holt) H,(tw,) do, des, (8.12) 


Lo 


E vt vt) vt) = — 
1 ? ã (27 )2 


Substituindo H (ww) pelas suas expressões obtém-se 


“4 
am) 


1 JO +64) 
E t) volt) vit) = A 2BRe) Flo std ETR Se 
VD vt) 7 E À Ao, to,) € | ado 


em que A é a área de qualquer dos hexágonos de KR, dada por 


à 
dA = 7 = (o + ow, —t,)? (8.14) 


Esta área é máxima para w. = w, + w,, atingindo então o valor 3 Aº/4. Nestas condições, teremos 


2 : midi 
E y,6t) y,(t) volt) = - Re | PAi), o t),) é hd | 
T a (8.15) 
expressão que nos vai permitir calcular W (mm mw). 
Com efeito, se fizermos 
8, 18,—84, = O (mod 25) (8.16) 
resulta 
g4* | 
E y(t) y,(t) yult) = Sa Re [Pulo + td) | (8.17) 
Se em vez disso fizermos 
Tr A—,=— E (mod 2) (8.18) 
vem 
34º 
E Ytt) Volt) Volt) Eae B-2 Im l Fm, : t3,) | (8.19) 


as 
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Conclui-se assim que a parte real e a parte imaginária do bispectro do sinal x(t), nas frequências 
ADA 09,) são dadas aproximadamente, à parte o factor 34º/8,º, pelos valores médios (que podem ser tempo- 
ra's, para os processos ergódicos) dos produtos dos sinais obtidos por filtragem estreita (ideal) de x(t) nas 
frequências (w,. (w, € w, + tw. € em que as fases (constantes) dos filtros obedeçam respectivamente às 
condições (8.16) e (8.18). (Por exemplo 4 =69,=6,=0, para a parte real, e 9 —- 4,=0, 9,=) /2, 
para a parte imaginária). 

Estas filtragens, análogas às estudadas para o espectro de potência pressupõem portanto filtros 
fisicamente irrealizáveis e como tal só são possíveis por via matemática, utilizando o computador. 

Notemos que a utilização de filtros realizáveis para a obtenção aproximada do bispectro por meios 
analógicos não se faz com a mesma facilidade com que se efectua para o espectro de potência, Por isso 


se preferem métodos aritméticos. 


9. COMPUTAÇÃO DO ESPECTRO DE POTÊNCIA E DO BISPECTRO 


Limitar-nos-emos a apresentar alguns tópicos dos métodos usados (1). 

Na prática, os métodos de filtragem numérica indicados no parágrafo anterior não são os mais indi- 
cados. Até 1965 eram os únicos utilizados para o bispectro (2). No que respeita o espectro de potência, 
p-rém, já antes disso se preferiam processos baseados em métodos estatísticos que partiam da estimação 
prévia da correlação, e que procuravam minimizar a variância das determinações (3). 

A partir de 1965, com a aparecimento do algoritmo — FFT — Fast Fourier Transform (4), foi possí- 
vel reduzir drasticamente o tempo de cálculo e realizar estimativas baseadas no cálculo do periodograma, 


dado por 


1 | 
| Chata == [2 Lua), [A (9.1) 
em que 
N—1 
= 
X(w)=4 > x(paye Po (9.2) 
p=o0 


sendo x(t) uma amostra de duração T= NA. A utilização apenas dos valores de x(t) intervalados de A 
implica a suposição de um espectro de banda limitada a | 3 | = q/à (teorema da amostragem). 

Este periodograma é depois ponderado por uma função J,(w) devidamente escolhida, denominada 
«janela» na linguagem da análise espectral, resultando finalmente para estimador do espectro de potência 
a expressão 


N— 
Es ns essi 8.3 
& (19) > n(u-TD) a (SÊ) (9.3) 
Ez 
tendo-se ainda 
aq/A 
Palo) | 
var É (ue | 32 (A) dá (9.4) 
2rT 
— q /A 


para espectros P(t)) variando lentamente dentro da largura da janela (5) ou seja, dentro da banda em que 
esta toma valores mais apreciáveis. 


Dado que J.(tw), além de ser par, apenas tem de obedecer à condição 


q/A 
a Int) dy =1 (9.5) 


o 


— m/A 


() Ver por exemplo, para o espectro, Jenkins & Watts, 1968; Brillinger, 1975; e para o bispectro, Brillinger & Rosen- 
blatt, 1967. 

(21 Ver Hasselman et al, 1963; MacDonald, 1965. 

Gi; Jenkins & Watts, 1968; Vários Autores, Technometrics, 3, 1961. 

(9) Para um estudo pormenorizado deste método, aconselha-se a ler Brigham, 1974. Quanto ao uso do periodograma 


já era conhecido antes de 195, sendo equivalente ao uso da correlação. Era porém mais lento que aquele, até à descoberta 
da FFT que lhe permitiu passar a ser o método mais rápido. | 
(5) 4 definição exacta de largura depende dos autores, mas anda À volta do conceito de concentração indicado. 
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a variância de (cw) pode reduzir-se arbitrariamente à custa do aumento da duração 'T (o que não acontece, 
prova-se, para o periodograma, que é uma estimativa não ponderada de P(r,)). Para T fixo, apenas se pode 
baixar a var (tw) à custa do aumento da duracão de J(wm)(!), o que vai prejudicar a resolução das 
estimativas, distanciando as frequências para as quais elas são indeperdentes. 

A redução do tempo de cálculo da expressão X(«w), através da FFT é espectacular, pois é da ordem 
de (1/N) log,N (por exemplo para N = 1000, a redução é da ordem de 1/100). 

Esta circunstância permitiu estender para o bispectro os métodos acima apresentados, com as adapta- 
ções necessárias. 


Assim, começa-se por calcular o periodograma de 3.º ordem 
(to, o» t4) =X (0) X WD) K(— w,—t,) (9.6) 


e daí, através das janelas bispectrais J.(w,,,), obedecendo à condição 


r/A 


1 
Dr ff Telmo, » t),) dety, dty, = 1 (9,7) 


. — m/A 


obtém-se para estimador do bispectro a expressão 


N-1 N-à 


A 1 2x5, 2 ms 2Ts, 2 TT S4 
+ (w, ' o) = 7 >, > Jal om T , L9g— T Em [a : To (9.8) 


s4—0 s=0 


tendo-se ainda 


1 


var Y (tw, td,) = (CRT 


FIA 
b (13,) P (19,) E (mo, + 00,) ff Ty (1,4) dá dy (9.9) 


para espectros variando lentamente dentro da largura (em cada dimensão) da janela. 

As considerações feitas para o espectro de potência a propósito do compromisso entre variância e 
resolução das estimativas põôem-se também para o bispectro. Pode provar-se que, em igualdade de reso- 
lução, o que se traduz por igual ordem de grandeza entre a largura (em cada dimensão) de J,(w,,w) ea 
largura de J (ww), e para uma série x(p4) normada, isto é, com variância unitária (2), se tem aproxima- 
damente, para o ruído branco (&P(w) constante) 


A 
de (1), + tds 
Ainda (9.10) 
var (e) 


em que M é o número de bandas independentes (em cada dimensão), ou seja, o cociente entre a banda do 
processo (q/4) e a largura da janela, 

Assim, por exemplo, se considerarmos um processo limitado a 1000 Hz tem-se como ordem de pgran- 
deza var WY/var P = 16, isto é, as flutuações das estimativas do bispectro são aproximadamente 4 vezes 
mais elevadas que as do espectro de potência. Quanto maior for o número de bandas, maior será este agra- 
vamento das variações das estimativas do bispectro em relação às do espectro de potência, 

Esta maior instabilidade resulta da existência de 3 factores aleatórios X (mw) em vez dos 2 usados 
para o espectro de potência, e constitui uma séria limitação das aplicações do bispectro, pois a redução da 


(1) A forma da função J,(w) também influi na variância, Para as janelas mais usuais, como as de Tukey, Farzen, Pa- 
poulis, a diferença não é porém grande. 
(?) A normação consegue-se substituindo a sucessão X, = X(DÃ) por x, / Va, sendo x,2 o valor quadrático médio 


calculado da amostra. Os valores x, supõem-se além disso préviamente deduzidos do seu valor médio, de modo a ter-se 


x = 0. 
a 
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variância exigirá um número muito elevado de dados. Este facto vai onerar o cálculo, pois muito embora a 
computação dos factores X(r)) seja extremamente facilitada pela FFT, a efectivação das médias ponderadas 
a 2 dimensões vai implicar um número de multiplicações que é sensivelmente da ordem de Nº, ao passo 
que para o espectro de potência aquele número é da ordem de N. 


Esta limitação não é, no entanto, o maior obstáculo, dada a crescente velocidade de cálculo dos 
computadores digitais. A maior limitação ao uso do bispectro reside na própria dificuldade em obter nos 
casos práticos um número suficientemente elevado de dados, em virtude da natureza limitada dos registos 
usualmente disponíveis. 


Não queremos terminar este artigo com uma nota pessimista. Na realidade hã que ter em conta a 
grande sensibilidade do bispectro na detecção dum processo não gaussiano. 


Para o ilustrar vamos apreseriar dois exemplos, através de gráficos obtidos na determinação do 
bispectro de processos simulados em computador. 


O primeiro (1) refere-se a uma situação em que estão presentes um processo gaussiano (ruido) e um 
processo constituido por uma sobreposição de impulsos de amplitude aleatória, ocorrendo também aleato- 
riamente no tempo segundo uma lei de Poisson (sinal). Apresentam-se duas determinações correspondentes aos 
casos só ruído e ruído + sinal, com a relação S/N = 20 GG. Conclui-se que, apesar de o ruído ser 5 vezes 
maior que o sinal (em potência), é visível a presen;a deste através do bispectro (Figs. 9.1 e 9.2). 


Teárico kr= — TEÓFICO Si 
sExperimental S/N 0 ” Experimental /N = 0.2 


Im p [to 635) 


Fig. 9.1 Fig. 9.º 


O segundo (2) refere-se a uma situação em que está presente um processo não gaussiano resultante 
da passagem de ruído-branco gaussiano estacionário por um sistema com uma não-linearidade de 2.º ordem, 
dado por y =x + qxº. O sinal de saída y(t) é normado de forma a ter-se y=0 e y2=1. Devido a proble- 
mas de amostragem resulta um fenómeno de sobreposição de bispectros que conduz a um valor teórico 
constante (2). Comparem-se (Figs. 9.3, 9.4) os casos q = O (caso linear) e q — 0,1. A pesar da pequenez da 
não-linearidade, esta torna-se também perfeitamente visível através do bispectro. 


e A 


(1) Ver Rebelo Simões, 1971. 
() Ver Rebelo Simões, 1972 e Rebelo Simões & Miranda Mourão, 1972, 
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e Experimental 


— T2ÓFico 
e Experimental 


Im Plcoçco) 
- | Ni . 


0, 
Q 
-1.5 
-1,0 
1234 %67BI7HHWMNR 
Fig. 9.8 Fig. 9.4 


Em qualquer dos casos os pontos calculados correspondem ao sector representado na Fig. 9.5, com 
N N 
tw), Sd O 2” tendo-se feito M = 8. 
MA * MA 


Os pontos foram ordenados arbitrariamente da forma indicada na figura e, para efeitos de represen- 
tação do bispectro, colocados como abcissas a intervalos iguais. 

Em conclusão, fica patente através destes exemplos a forma notável como o afastamento da distri- 
buição dum processo em relação à gaussiana transparece na determinação do bispectro, mau grado as 
flutuações das estimativas. 
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APÊNDICE 
TRANSFORMAÇÃO LINEAR DUMA DISTRIBUIÇÃO MULTINORMAL 


Consideremos as variáveis (X,,X,,-.. X,) obedecendo a uma distribuição conjunta gaussiana ou 
multinormal 


- NON 
1 : 
piso aa ad E E 
ja 9 


N;/2 1 
(297) ; (det a) nl  l=l 


VA) (x—m.) | (A1l)=(Tl) 


Vamos usar, para facilitar a análise, a notação matricial 
Assim, seja X a matriz coluna das variáveis x, M a das médias m,, e façamos 


a T=X—M (A 2) 


o que equivale a definir as novas variáveis y =x —m,, de média nula. 
A expressão (A 1) da distribuição multinormal tomará a forma abreviada (1) 


pt pm = enaga | rã (A 3) 


N/2 1/2 
(25) (det A) 


ma 


Por outro lado a expressão (7.3) À, = E(x,—m,) (x —m,) pode escrever-se 


AS EYY (A 4) 


Consideremos agora uma transformação linear das variáveis y, definida por uma matriz quadrada A de 
ordem N 


Y=AY (A 5) 
equivalente a 
N 
Ro x a d=L8 aan (cf 7,4) 
Ter-se-á então (3) ee 
EY' = ABY = 0 (A 6) 


pelo que as correlações de 2.º ordem das variáveis que definem Y coincidem com os momentos, constituindo 
a matriz A' dada por 


E T 
N= EYT'=EAT(AT) = ANA. (AT) 


em virtude de (Ad e A 5). 
Esta relação mostra como se obtêm as correlações das novas variáveis y', em função das correlações 
das antigas y,. 
Resta determinar a distribuição das novas variáveis y'.. 
Para isso basta exprimir o argumento da exponencial em (A 3) nas novas variáveis, e atender à 
identidade 


(4) 
PO Ye Y'y) dy”, dy, dy =P (07,00 Yy) dy, dy, .... dyy (A 8) 


pelo que se tem, como se sabe do cálculo, 


A ; o(y pi gera: De) 
PO o e D) =P og ceB) | O | (A 9) 
Io VD | 


T 
(1) A =matriz transposta de A 
(*) Seria mais correcto escrever p,(x) = p,ly), com T = X-M, mas omitimos, para facilitar a notação, os indices, que 
se subentendem. 
0) O zero representa aqui a matriz coluna constituída por zeros. 


(9) A transformação linear (A 5) implica uma solução única Y = AE para o sistema de equações correspondente, 
pelo que (A 8) só tem um termo. 
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em que 9(Y,,Y,--Yy)/9(9,,7',:-- Y'y) é o Jacobiano da transformação (A 5), dado pelo determinante 
da matriz J, de elementos 


= à Yi (A 10) 
dy" 
Ora de (A 5) tira-se 
-— ] 
TesA TWT (A 11) 
se a matriz A não for singular 
Da transformação linear (A 11) obtém-se 
O 
Yi = 8", (A 12) 
os”, 
é —1 
em que a”. são os elementos da matriz A. 
7 a. | 
Daqui se conclui que J = À (A 13) 
«id | 
donde det J= (det A) (A 14) 
pelo que (A 9) se pode escrever, abreviadamente 
p(T) =p(Y)/|det A | (A 15) 
T 
Calculando agora o produto Y ['Y, argumento da exponencial (A 3), vem 
T -1 T=-=1 =1 T 2“T-=1 =1 inca T =1 7 
TIT=(i Ti A T=Tr (A! x AYT=rimMaA) EW=TPT (A 16) 
em que p= (a) (A 17) 
Por outro lado, de (A 7) tira-se 
det na'=det A (det A)2 (A 18) 
donde finalmente, atendendo a (A 3), (A 15), (A 18), vem 
(77) ; À 
À TD e ss ones SUA ' ' 
Pp N 12 p 3 To pr Y (A 19) 


(2a) a (det a') 


que mostra que as novas variáveis y' resultantes da transformação linear (A 5) possuem igualmente uma 
distribuição multinormal, sendo a matriz das correlações dada por (AT). 

Este resultado é muito importante, pois permite-nos, por uma escolha conveniente da transformação A, 
simplificar notavelmente a expressão de distribuição. 

Assim, é possível obter um novo conjunto de variáveis cuja matriz de correlações seja diagonal. Na 
realidade, sendo A simétrica, basta usar na transformação Y = AX uma matriz A ortogonal (1) cujas 
linhas sejam constituídas pelos vectores próprios normados de a. 

Nestas condições, a nova matriz a' será diagonal, sendo os seus elementos não nulos os valores 
próprios de a. 


() Se os vectores linhas são ortogonais entre si e normados, tem-se Z 4,, 8, = 
Es 


T 
Fazendo a,,=b,,. tem-se E a; by=3,,» ou seja AB=I. Mas B=A . 
k 


| T 
Donde B=A =A que define a matriz ortogonal. Fácil é de ver que também os vectores colunas resultam orto- 
normados. 
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Com efeito, sendo À e X, um valor próprio e um vector próprio de a, tem-se 
AR = A, (A 20) 


Por outro lado prova-se sem dificuldade que, sendo A simétrica, os seus valores próprios são reais, 
e a valores próprios distintos correspondem vectores próprios ortogonais (1). 


Assim, sendo 
am a (A 21) 


a matriz A calcula-se facilmente. 


E = - or o ” 
A=A(A AY = x! DD ese (A 22) 
ii | | (O Ag m=e======= () 
! ] | WU 7 
| | | 
| E E ày Xy = 4! 
T | | 
Ay OQ esses: Xe 


dada a ortonormalidade dos Ke 


Ou seja 
Nx = Md (A 28) 


o que prova que a matriz de correlação A é de facto diagonal. 
Se introduzirmos agora este resultado (A 22, A 23) na expressão geral (A 19) da distribuição p (Y'), e 
atendermos a (A 17), obteremos a forma simplificada pretendida. 


Assim, 
1 1 ' a 
PV) =P(y Se y)=—— .———— > exp A > di (A 24) 
N/2 1/2 E) " 
(257) (A goes Dad 5 fi El 
que se pode escrever sob a forma dum produto de N factores semelhantes 
N 
PSF Yp= HH PY) (A 25) 
E-1 
em que 
(y') a fimo uns DO 
BR = e: (RS fe re 
i (27A,)1/2 p 2 À (A 26) 
Notemos finalmente que o facto de A' ser uma matriz de correlações nos mostra que se tem 
M=hasEy;i>o (A 27) 
Os valores próprios de 4, além de reais são, portanto, positivos. 
Por outro lado, do carácter ortogonal de A resulta 
det A=det A=A, À,...Ay>DO (A 28) 


Assim, o determinante da matriz das correlações mantém-se invariante sob uma transformação 
ortogonal e é positivo, 


(1) Se se der o caso de n valores próprios serem iguais, ao valor comum corresponderá um espaço vectorial de di- 
mensão n, donde é possível obter n vectores ortogonais entre ai. 
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Em homenagem ao Professor Aureliano de Mira Fernandes 


Sobre a elasticidade relativista dos corpos rígidos 


| — À lei de Hooke relativista 


RESUMO 


Prepõe-se a introdução nos cursos de iniciação 
à Teoria dz Relatividade das fórmulas da elasti- 
cidade relativista dos corpos rígidos, estabelecidas 
por Mec Crea por processos elementares. 

Os paradoxos relacionados com o movimento 
dos corpos rígidos em Relatividade (paradoxo do 
disco a rodar e paradoxo dos 3 grous de liberdade 
dos corpos rigidos em Relatividade Restritiva) podem 
ser esclarecidos desde que se faça aq distinção entre 
os conceitos de rigido e indeformável, 

É apresentada a equação do movimento (com 
deformação) de uma barra rígida não homogénea. 


A «elasticidade relativista dos corpos rígidos» é uma matéria 


ANTÔNIO BROTAS 
Instituto Superior Técnico 
Universidade Técnica de Lisboa 


SUMMARY 


We propose the introduchon in lhe elementary 
courses about The Theory of Relativity of the rela- 
tivístic formulas of rigid bodies elasticity esta- 
blished by Me Crea by elementary processes, 

The paradoxes relationed with the mouvement 
of the rigid bodies in the Theory of Relativily (lhe 
peradox of the sliding disc and the paradox of the 
three degrees of freedom of rigid bodies in the 
Special Theory of KRelativity) can become clear df 
we disting the concepts of rigid and indejformable. 

We present the equation of the mouvement 
(with deformation) of a non-homogeneous bar. 


de todo ignorada nos tratados de 


Relatividade (*) e só referida numa meia dúzia de quase desconhecidos artigos. A sua inclusão nos cursos 
de iniciação é, no entanto, possível e inteiramente natural. Hstou certo que, quando se generalizar, fará 
desaparecer pela base um conjunto de «paradoxos» actualmente considerados como ainda não esclarecidos 
pela maioria dos relativistas. 

A nota que se segue é facilmente compreensível para quem conheça a Teoria da Relatividade, tal como 
há uns anos atrás era ensinada nas cadeiras de Física e de Mecânica Racional do I. 8. T. Ao apresentá-la, 


faço votos para que, a breve prazo, surja na escola o interesse suficiente para que este tipo de matérias 
possa, de novo, nela ser ensinado, 


Vamos admitir que corpo rígido é aquele em que 
as ondas de choque se propagam com a velocidade 
c (ver nota 1, página [771]. 


Consideremos uma barra rígida, homogénea, de 1) 


S=S Pl=hVi=p (*- = 


Cc 


secção S,, comprimento |, e densidade d no refe- > ê aa E uá 
rencial próprio. E é Mo É =" Sa la Eq a Y Sa do Po 
Vamos supor que esta barra se desloca num E E É ==: ER 


referencial A (ver nota 2, pág. [77] com uma veloci- 
dades v. Neste referencial A a secção, o comprimento, 
a energia e a quantidade de movimento da barra são 
dados, respectivamente, por: 


Vamos admitir que num dado instante t (t — tem- 
po de referencial A) uma extremidade da barra, a 
extremidade da frente, por exemplo, pára devido ao 


(*) €C, MOLLER, por exemplo, no seu tratado «The Theory of Relativity> dedica algumas poucas páginas ao estudo 
dos sólidos elásticos. À pág. 179 (da edição de 1962) escreve: «O tensor das tensões está relacionado com a deformação 
interna da matéria. No referencial próprio esta conexão é dada pelas equações da teoria da elasticidade usual; assim, 
para as pequenas deformações é dada pela lei de Hooke>s. Falta, pois, o essencial, a determinação da lei elástica relativista. 

Manuscrito recebido para publicação em 7/17/78. 
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choque com uma parede indeformável, ou com uma 
barra igual vinda em sentido contrário com a mesma 
velocidade v (ver nota 3 página [77]. A segunda 
extremidade da barra continua em movimento com a 
velocidade v. Uma vez que supomos o material rigido, 
o efeito do choque propaga-se com a velocidade c. 
Um cálculo elementar mostra-nos que a extremidade 
que continuou em movimento só é atingida por este 
efeito no instante 


VIP 


|] 


ELA] Ig = E, n 


C W 


Vamos admitir que a onda de choque, que se 
propaga ao longo da barra, separa uma parte que 
fica parada e comprimida da parte não comprimida e 
ainda em movimento. No instante t,, a barra está 
inteiramente parada (em A) e o seu comprimento é 


4) Le di l, V 1 = pf : cá : 
Cc+-v 


Se admitirmos que a barra é elástica o choque 
deve ser reversível. Neste caso a partir do instante t,, 
a segunda extremidade deve voltar para trás com a 


velocidade v, e no interior da barra uma onda de uno- 
que (de sentido contrário à anterior) deve separar a 
parte não comprimida, e de novo em movimento, da 
parte comprimida e em repouso, 

No instante t,=2 t,-t, a onda de choque atinge 
a extremidade que inicialmente foi parada e a barra 
volta, de novo, a estar toda em movimento com a 
velocidade v, mas agora em sentido contrário. 

Durante o intervalo [t, t,] a extremidade da 
barra que chocou mantém-se encostada à parede (ou 
à outra barra) que sobre ela exerce uma força que 
representamos por F. 

sendo o choque elástico e reversível a energia 
mantém-se constante, 

Concluímos, assim, que a energia da barra parada 
e comprimida com o comprimento 1, deve ser igual 


2 
po 


V1i—s 


com a velocidade v. 


à energia E — da barra não comprimida e 


Esta informação é suficiente para determinar a 
lei elástica do material rigido. 

Imaginemos, com efeito, um processo lento de 
compressão de uma barra (imóvel em A) desde o 
comprimento 1 até ao comprimento 1. 


TS E TM 

v 1 I 
eim] 
| I 

Di J 


Fig. 1 


(4 
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ul + E” 
| CO. 
cu pda : 
º| e 
Le a? 
1d 
E 
o [o] 
- E 
MI 
h > 
Fig. E 
O trabalho dispendido deve ser igual ao aumento Fazendo 
d: energia 
le ! c 
| É 4 5) S, = -Vi-s E 
4) Es E = me mta pas À lo Cc-+ vw 
vi—ya 
| 
I 
| 
| P 
I 
I 
l 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
l 
| 
| 
| 
| 
I 
I 
| q 
! compressão N q a tracção 
s=0 551 E 5 —» 09 


— oo e ae e. o o ic Cc = ic — cc — cs O o EO CE E c—i—Ooc— — o —o— —O — —O 


Pig. & 


TÉCNICA 449/450 To 


e usando as relações 


6) É EE === Ef 


podemos escrever 


2 
Mm, € l à 
Se Ê = — 
p Se 


Por outro lado, sendo p a pressão na secção 
transversal da barra, variável durante a compreensão, 
podemos escrever 


7) z = 


Tendo em conta que as duas últimas expressões 
de = devem ser iguais qualquer que seja o valor de 


8.» derivando, obtemos 


E (= a) —Sbr(so 


o “dl 
o a &E(L-) 


E esta a lei elástica dos corpos rígidos, lei não 
linear, que podemos considerar como a correspondên- 
cia relativista no caso limite dos corpos rígidos da 
lei de Hooke (ver nota 4, página [77]. 

Como seria de esperar, do lado da compressão, p 
tende para infinito quando s tende para zero, 

Do lado da tracção é curioso verificar que p 


ad ed 


tendo assintoticamente para 


quando s tende 
para infinito (ver nota 5, página [77]. 

Nesta concepção, os corpos rigidos aparecem, 
assim, como corpos que se podem esticar indefini- 
damentea sem nunca atingir uma tensão infinita. 


Representamos por “Da a densidade de um corpo 
rígido, no referencial próprio, quando não compri- 
mido. Tem interesse conhecer 6, — 5 (5), densidade 


do mesmo corpo, no referencial próprio, quando 
deformado. 
A relação 
- Mg e 
11) Ê =. ———+*—— ps EE cg 
V1i—p 
76 


permite-nos, facilmente, chegar ao resultado 


12) “h= bo ( E E) 
2 sº 


Quando s tende para infinito, p, tende para 
um valor que é metade da densidade inicial. 


Indiquemos agora um segundo caminho para 
chegar ao resultado 10), que pode ser olhado como 
uma confirmação da hipótese que fizemos sobre a 
propagação das ondas de choque no interior da 
barra. 

Consideremos um intervalo de tempo àt contido 
no intervalo [t,,t,]. Representemos por Am, a quan- 
tidade de massa própria (não comprimida) que no 
intrevalo A t passa do estado de movimento ao repouso 
tem A). 

Um estudo cinemático relativamente simples per- 
mite-nos escrever 


Ss(v+Hev tm 
V1-8 


A variação da quantidade de movimento da barra 
neste intervalo de tempo At é dada por 


13) ATA) = 


14 4 am, V Siv+te)v PAL 
+] É o E 
| Vip 1— 

dP 


, em que f é a resultante de 


A equação f= E 
todas as forças exteriores que se exercem sobre um 


dP 
corpo e. à -a sua variação de quantidade de movi- 
mento, é válida em Relatividade tal como em Física 
Clássica. 

Durante o intervalo de tempo At a única força 
que se exerce sobre a barra é a força que atrás 
designamos por F. Temos ainda, manifestamente, 
F=p.S, sendo p a pressão nas secções da barra já 
comprimidas. 


Reunindo estes elementos podemos escrever 


1 dP 
S dt 


te (1.4) 
2 sº á 


o que condiz com o resultado anterior. 


viv+o) 


É AO 6 as 
| — pº 


Ra 


15) p= 
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NOTA TERMINAL 


Os resultados expostos, embora muito simples e 
baseados em noções elementares de Relatividade, só 
foram publicados pela primeira vez em 1952, nos 
artigos de Mc Crea citados na bibliografia. No pri- 
meiro artigo, Mc Crea deduz a equação da relativi- 
dade dos corpos rígidos a partir da conservação da 
energia e no segundo, em colaboração com Hogarth, 
faz uma verificação usando a conservação da quanti- 
dade de movimento, 

Tive conhecimento destes artigos por intermédio 
do artigo de Cavalleri, numa altura em que estudava 
o disco a rodar. Convencera-me de que o problema do 
disco era um problema de elasticidade, mas a duas 
dimensões. Tinha chegado aos mesmos resultados de 
Mc Crea (só que encontrados por ordem inversa, pri- 
meiro a partir da quantidade de movimento e depois 
da energia) e tinha encontrado uma generalização 
para o caso a duas dimensões que não considero intei- 
ramente satisfatória. Na ocasião escrevi a Mc Crea a 
perguntar-lhe se tinha avançado para além dos arti- 
gos de 52 e se tinha utilizado as suas fórmulas no 
disco a rodar. 

Respondeu-me, com grande gentileza, dizendo 
que não, mas que, como eu bem suspeitara, sempre 
tinha continuado a pensar no problema do disco a 
rodar. 

Este assunto pode, eventualmente, ter interes- 
santes possibilidades de desenvolvimento. Um deles 
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(que se afigura muito duro) é o estudo das ondas 
transversais nos corpos rígidos. O estabelecimento 
de convincentes fórmulas de elasticidade a 2 e 3 di- 
mensões terá ,eventualmente, de ser feito em conju- 
gação com o estudo destas ondas transversais. 

No caso a uma dimensão, a equação que julgo 
inédita 


ap EN see 
16) x|- AE Elo, +1] 
1 x O X 


(em que x e t são as coordenadas do referencial de 


inércia e x a coordenada associada aos diferentes 
pontos duma barra e correspondente às distâncias 
medidas sobre a barra quando não deformada) é a 
equação do movimento longitudinal de uma barra 


v (x) variável. 


O estabelecimento desta equação exige bastante 
trabalho e muito cuidado no manipular das compo- 
nentes do tensor impulsão energia expresso nos 
referências de inércia correspondentes em cada ina- 
tante a cada ponto da barra, Mas exige, sobretudo, 
a convicção inicial de que se vai chegar ao mesmo 
resultado quer se utilize a conservação da energia 
quer se utilize a conservação da quantidade de mo- 
vimento. 


rígida não homogénea: 9! = 


[4] A. Brotas — C. R, Acad. Sc. Paris t 267 (1968). 
[5] A. Brotas — L. Nuovo Cimento (Fev. 1969). 
[6| A. Brotas — These, Paris. A. O. 3081. 


1 O corpo rigido relativista é, portanto, nesta concepção um corpo deformável. Foi, exactamente, a identificação 


entre rigido e indeformável que, trazida da Física Clássica para a Relatividade, deu origem a uma série de paradoxos 
relacionados com a dificuldade em compreender o movimento dos corpos rigidos (e por extensão dos sólidos em geral) 
em Relatividade. O mais célebre destes paradoxos é o paradoxo do disco a rodar segundo o qual em Relatividade Res- 
trita um disco indeformável que esteja a rodar não pode parar e um disco que esteja parado não pode entrar em rota- 
cão. Estas dificuldades são aspectos particulares do paradoxo mais geral, posto em evidência por Born, que consiste no 
facto de em Relatividade Restrita os corpos rigidos-indeformáveis terem só 3 graus de liberdade. Em face desta difi- 
culdade a quase totalidade dos autores que abordam o assunto abandonam o estudo do problema no âmbito da Relati- 
vidade Restrita e procuram encontrar no quadro da Relatividade Generalizada, à custa das deformações do espaço- 
“tempo, os 6 graus de liberdade que não encontram em Relatividade Restrita. Os trabalhos deste tipo tornam-se, em geral, 
extremamente confusos e parece não terem conduzido a nenhumas conclusões claras. A via preconizada neste artigo é a 
oposta, A Identificação entre indeformável e rígido é quebrada. Indeformável é olhado como uma noção geométrica. 
Rígido como uma nocão física. Corpo rigido significa simplesmente o corpo mais rigido possível. Em Fisica Clássica, 
nada se opõe a que consideremos ser indeformável o corpo mais rigido possível. Em Relatividade, temos a imposição 
das ondas não se propagarem com velocidade superior à da luz. Daqui a necessidade do corpo rigido ser deformável. 

2) Sempre que neste texto falamoz em referencial deverá ser entendido referencial de inércia. A noção de refe- 
renciais não de inércia da Física Clássica estava associada à noção de corpos rigidos que embora sujeitas a forças ne 
não deformavam. Como em Relatividade excluimos a ideia de corpos rigidos-indeformáveis temos também de por de 
lado a noção de referenciais não de inércia. 

3) A barra a vir do outro lado é um artifício para evitar falar em parede indeformável, 

Neste texto escrevemos, em geral, compressão mas todas as fórmulas e raciocínios são válidas para as tracções. 
Há, porém, que imaginar um dispositivo mais sofisticado para parar a extremidade de traz da barra. 

4) No caso de um material não rígido limite as ondas de choque devem-se propagar no material não comprimido 
e no referencial próprio com uma velocidade a inferior a c. A partir deste dado pode ser determinada a lei elástica. 


db) No limite, em valor absoluto, a relação PD. é igual à relação entre a pressão e a densidade de uma onda elec- 
| 
Ps 
tromagnética. Neste caso o sinal vem, porém, trocado porque p é uma tracção. 
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CURSOS DE FORMAÇÃO EM HIDROLOGIA 


A — INTRODUÇÃO 


O C. E. S.U.R. (Centro de Sistemas Urbanos e Regionais das Universidades de Lisboa), os departa- 
mentos de Hidráulica e de Urbanização e Sistemas do I.5.T. julgaram oportuno realizar em colaboração 
com a Direcção Geral dos Recursos e Aproveitamentos Hidráulicos dois Cursos sobre Hidrologia. Estes cur- 
sos não têm características de curso intensivo e têm como objectivo contribuir para uma preparação cien- 
tífica de base e conhecimentos técnicos actualizados que permitam aos licenciados interessados em Hidro- 
logia — no domínio do palneamento, do projecto ou da gestão — utilizar convenientemente as potencialidades 
quer da Hidrologia Física, quer da Hidrologia Estatística. 

Durante estes cursos serão apresentados casos de aplicação e os participantes terão a possibilidade 
de utilizar o Centro de Cálculo das Universidade de Lisboa e as bibliotecas do C. E, S. U. R. e do 1. 8. T. 


B — DESCRIÇÃO DOS CURSOS 


CURSO SOBRE HIDROLOGIA FÍSICA (HF) 


1. PROGRAMA: 


— Descrição física dos processos hidrológicos 
— Variáveis hidrológicas 

— Medições hidrológicas e redes de observação 
— Processamento, crítica e publicação de dados 
— Modelos determinísticos 


2. DOCENTES 


Director: António Quintela (I. S.'T.); Pinto Peixoto (Fac. de Ciências); Zózimo Rego (TI. 5. A.); Emi- 
dio Santos (I.S.T.); João Hipólito (L.S.T.). 


CURSO SOBRE HIDROLOGIA ESTATÍSTICA (HE) 


1. PROGRAMA: 


— À teoria do aleatório e os tipos de variabilidade hidrológica 

— Fundamentos da Teoria das Probabilidades e da Informação úteis à Hidrologia 
— Descrição, filtragem e síntese estatísticas em Hidrologia 

— Inferência e Amostragem 


2. DOCENTES 


Director: Luís Valadares Tavares (I. 8. T. e C E.S.U.R); Técnicos da D.G.R.A.H.: F. Nunes 
Correia (I.S.T.); J. E. Silva (1. S.T.); A. Porta Nova (C.E.S.U.R. e I. S.T.); Delmar Baptista 
(CE.S.U.R. e LS. T.) 


3. CONDIÇÕES GERAIS: 


Estes cursos realizam-se no C.E.S.U.R. (I.S.'T., Av. Rovisco Pais, Lisboa) com o seguinte horário: 
9h30 - 12h30. 

O curso HF decorre de 2-11-78 a 17-11-78 e o de HE de 20-11-78 a 6-12-78. 

Quaisquer esclarecimentos adicionais podem ser obtidos através de: L, Valadares Tavares (C. E. S. U.R,, 
Instituto Superior Técnico, Av. KRovisco Pais, Lisboa. Tel. 882992). 

Os partic pantes que desejem obter um certificado de aproveitamento podem requerer, para tal efeito, 
a realização de uma prova que terá lugar a seguir ao final da cada curso, 
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Resumo dos Artigos publicados na «Técnica» n.º 449/450 
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C. D. U. 51.851 


MONTEIRO, Antônio Aniceto — Conjuntos gradua- 
dos de Zadeh, «Técnicas, Lisboa, LIIl (449/450) 
Set. /Out. 1978, p. 11-34. 


C. D. U. 538.3 


FARO, M. Abreu — Sobre o estabelecimento das 
condições fronteiras de natureza electrodinâmica 
numa superficie em movimento, «Técnica», Lisboa, 
LI (449/450) Set./Out. 1978, p. 35-44, 


É demonstrado que as condições fronteiras numa 
superfície em movimento se podem obter a partir das 
condições fronteiras dos corpos em repouso utilizando 
a Transformação de Lorentz que imobiliza num deter- 
minado referencial o ponto em movimento da referida 
superfície. 

Para tanto, foi necessário o prévio estabelecimento 
das leis de transformação das densidades de carga e de 
corrente superficiais, o que foi feito com generalidade. 

As condições fronteiras estabelecidas coincidem com 
as condições fronteiras obtidas a partir da forma inte- 
gral das equações de Maxwell e utilizando os teoremas 
de Helmholtz associados a superfícies e contornos de- 
formáveis no tempo, 

Na demonstração postula-se que determinados inte- 
grais tendem para zero o que é fisicamente aceitável 
mas obriga a analisar com cuidado os modelos em causa 
e as leis de variações expectáveis das densidades de 
carga e de corrente, 


TÉCNICA 449/450 


Ano LIll — Setembro/Outubro 1978 


+ 


C. D. U. 311.15:535.33 


SIMÕES, F. E. Rebelo — Alguns aspectos da corre- 
lação de 3.º ordem e do bispectro dos processos 
estacionários, «Técnicas, Lisboa, LII (449/450) 
Set./Out. 1978, p. 45-71. 


Introduzem-se, como generalização da 2.º ordem, os 
conceitos de correlação de 3.º ordem e de espectro de 
3.º ordem (bispectro) dos processos estacionários, estu- 
dando-se as suas simetrias. Faz-se em seguida idêntica 
extensão para correlações e espectros cruzados. Esta- 
belecem-se as relações entrada-saída de sistemas linea- 
res invariantes para correlações e espectros, Considera- 
-se o caso particular dos processos gaussianos estacio- 
nários de média nula, para os quais se verifica que o 
bispectro é nulo, deduzindo-se daí as aplicações práti- 
cas desta análise, Interpreta-se a seguir o bispectro 
como resultado de filtragens lineares, à semelhança do 
que se faz para o espectro de potência. Finalmente in- 
dicam-se alguns tópicos sobre a computação do espectro 
de potência e do bispeciro, comentando-se os prós e os 
contras da análise espectral de 3.º ordem, 


C. D. U. 533.2:530.12 


BROTAS, António — Sobre a elasticidade relativis- 
ta dos corpos rigidos: a lei de Hooke relativista, 
«Técnica», Lisboa, LIII (449/450) Set./Out. 1978, 
p. 73-77. 


Propõe-se a introdução nos cursos de iniciação à 
Teoria da Relatividade das fórmulas da elasticidade 
relativiata dos corpos rígidos, estabelecidas por Mc Crea 
por processos elementares. 

Os paradoxos relacionados com o movimento dos 
corpos rigidos em Relatividade (paradoxo do disco a 
rodar e paradoxo dos 3 graus de liberdade dos corpos 
rígidos em Relatividade Restritiva) podem ser esclare- 
cidos desde que se faca a distinção entre os conceitos 
de rigido e indeformável, 

É apresentada a equação do movimento (com de- 
formação) de uma barra rigida não homogênea, 
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Synopsis of articles published in «Técnica» n.º 449/450 
Lil Seprtember/October 1978 


U. D. C. 311.15:535.33 


SIMÕES, F. E. Rebelo — Some aspects of the 3" 
order correlation and of the bispectrum of stationary 
processes, «Técnica», Lisboa, LIII (449/450) Set./ 
/Out. 1978, p. 45-71. 


As generalization of 2" order analysis, the concepts 
of 3rd order correlation and 3rd order spectrum (bis- 
pectrum) of stationary processes are introduced, and 
their symmetry properties are studied. Identical exten- 
sion for cross correlations and cross spectra is made, 
The input-output relations for correlations and spectra 
are established in linear time-invarlant systems. The 
gaussian case is considered in which the bispectrum is 
shown to be zero. Some practical applications of this 
fact are presented. Next the bispectrum is interpreted 
as a result of linear filtering operations, like power 
spectrum. Finally some topics on the computation of 
power spectrum and bispectrum are indicated, with 
comments on the advantages and limitations of 3rd 
order spectral analysis. 


U. D. C. 533.2:530.12 


BROTAS, António —- On the relativistic elasticity 
of the rigid bodies: relativistic Hooke's law, «Téc- 
nica», Lisboa, LIIl (449/450) Set./Out. 1978, 
p. 73-77. 


We propose the introduction in the elementary 
courses about The Theory of Relativity of the relati- 
vistic formulas of rigid bodies elasticity established 
by Mec Crea by elementary processes, 

The paradoxes relationed with the mouvement of 
the rigid bodies in the Theory of Relativity (the para- 
dox of the sliding disc and the paradox of the three 
degrees of freedom of rigid bodies in the Special Theory 
of Relativity) com become clear if we disting the 
concepts of rigid and indeformable. 

We present the equation of the mouvement (with 
deformation) of a non-homogeneous bar. 


o E O O 


— a 


U. D. C. 51.951 


MONTEIRO, António Aniceto — On [fuzzy sets, 
«Técnica», Lisboa, LIII (449/450) Set./Out. 1978, 
p. 11-34. 


U. D. C. 538.3 


FARO, M. Abreu — On the establishment of elec- 
trodynamic boundary conditions at a moving sur- 
face, «Técnica», Lisboa, LIIl (449/450) Set./Out. 
1978, p. 35-44, 


It is shown that the boundary conditions at a 
moving surface can be obtained from the boundary 
conditions of bodies at rest using the Lorentz Trans- 
formations which immobilizes on a certaln frame the 
point moving at the surface referred to. 

This required the previous establishment of the 
transformation laws of densities of surface charge and 
current, which was generally accomplished. The boun- 
dary conditions established coincide with the boundary 
conditions obtained from the integral form of Max- 
well's equations and using the Helmholtz theorems 
ERASeEADEM with surfaces and contours deforming with 
time. 

In all of this it was necessary to postulate thal 
certain integrals vanish, which is physically acceptable 
but requires a careful study of the models in question 
and of the which may by expected from the charge 
and current densities, 
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445-446  12-1-977/78 349 


62d — ENGENHARIA CIVIL 


Principios para a reestru- 
turação do ensino da 
Engenharia Civil ... ... 


C. D. U. 624:371.214 


445-446 12-1-977/78 173 


António de Caravlho Quintela 


O ensino da Engenharia 
Civil no nosso Pais . 
J. 8. Brazão Farinha 


Método directo de distri- 
buição de momentos 
para o cáleulo de estru- 
turas sujeitas a esforços 
horizontais ... ... 


Alfredo de Pinho Morgado 


Industrialização e quali- 
dade da construção ... 
José A. Teixeira Trigo 


C. D. U. 624:371.214 


445-446  12-1-977/78 169 


C. D. U. 624.04 


443 10-977 47 


C. D. U. 624.057 


4453-446  12-1-977/78 221 


C. D. U. 6240723341 
Vigas Vierendeel simples 447 2-978 397 
Jd. S. Brazão Farinha 


625 — CAMINHOS DE FERRO E ESTRADAS 


C. D. U, 625.8(661):691.223 
Les chaussées en sables 
gvpseux du Sahara ... 447 
José Carlos de O. 8. Horta 


2-978 417 


627 — CURSOS DE AGUA E CANAIS NATURAIS. 
RIOS, PORTOS, OBRAS PORTUÁRIAS E 
COSTEIRAS: BARRAGENS, ETC. 


C. D. U. 627.8.09(6/9.4) 

A engenharia nacional e 

a instalação de unida- 

des de grande potência. 

O caso de Cabora-Bassa  445-446  12-1-977/78 255 
Alberto Abecassis Manzanares 
Antônio Carrisso 
Antônio de Carvalho Quintela 
Antônio José da Costa Santos 
Jodo de Salvador Fernandes 


C. D. U. 627.82:627.89 
Desastres de barragens 
em 1976. O caso de 
Pe ton suis cerca via 443 10-977 7 
J. Laginha Serafim 
C. D. U. 627.89.007.1:174 
Acidentes em obras de 
engenharia civile a res- 
ponsabilidade dos enge- 
mheiros. O caso da bar- 


ragem de Teton ... ... 448 3-978 485 
Antônio José Correia Mineiro 
66 — INDÚSTRIAS QUÍMICAS, TECNOLOGIA 


QUIMICA 
665 — ÓLEOS. GORDURAS. CERAS 


C. D. U. 665.51/58 
O complexo petroquímico 
de Sines; concepção e 
estado actual dos tra- 
balhos ; 445-446  12-1-977/78 305 
Jodo Pedro Antunes Gomes 


TÉCNICA 449/450 


Artigo N.º Data Pão. 


C. D. U. 665.51/58 
Petroquímica ... .. 445-446 12-1-977/78 297 


Luis L. Ramalhete 


666 — VIDRO. INDÚSTRIAS CERÂMICAS, ETC. 
C. D-. U. 666.3/7 
Materiais cerâmicos: mi- 
croestrutura e proprie- 
CMOS quai mia dvd mes: ai 444 
C. A. Pacheco da Silva 
MH. Helena Carvalho 
M. Manuela Oliveira 


ll-BTTo TT 


669 — METALURGIA 


C. D. U. 669.1:543.42.64 
Determinação de coefi- 
cientes de repartição 
sólido /líquido em ligas 
ternárias Fn-Mn-C 445-446  12-1-077/78 355 
M. Helena de Melo Carvalho 


C. D. U. 669.1:621.762 
Pulverometalurgia de li- 
gas ferrosas ... ... ... 445-446 12-1-977/78 327 
Carlos Alberto Pacheco da Silva 


C. D. U. 669.14:621.785:620.163 
O controle de qualidade 
nos tratamentos tér- 
FOCO esses + 445-446  12-1-977/78 313 
Antera Faleriana de Seabra 


C. D. U. 669.14:621.785:001.2 
Necessidade de investiga- 
cão e de desenvolvi- 
mento tecnológico em 
Portugal na área dos 
tratamentos térmicos 
dos aços 445-446  12-1-977/78 345 


Antônio de Pádua Loureiro 


68 — OFÍCIOS, ARTES E INDUSTRIAS ESPECIA- 
LIZADAS PARA ARTIGOS ACABADOS E 
VESTUÁRIO 


681 — MECANISMOS, APARELHOS E MAQUINAS 
DE PRECISÃO. INSTRUMENTOS 


C. D. U. 681.14.016 
Uma bibliografia de 
APL — 2.º parte ... ... 447 


f. Marques Henriques 


2-978 431 


TÉCNICA 449/450 


Artigo N.º Data Pág. 


69 — CONSTRUÇÃO DE EDIFÍCIOS, MATERIAIS 
E PROCESSOS DE CONSTRUÇÃO 


C. D. U. 69.002:72,011.1 


Tecnologias da constru- 
ção de habitação ... ... 448 
José A. Teixeira Trigo 


S-D7B 4653 


691 — MATERIAIS DE CONSTRUÇÃO 


C. D. U. 691.223:625.8(661) 


Les chaussées en sables 
gypseux du Sahara ... 447 
José Carlos de O, 8. Horta 


2-978 417 


697 — INSTALAÇÕES DE AQUECIMENTO, VENTI- 
LAÇÃO E AR CONDICIONADO 


C. D. U. 697.7:621.47:517 


Métodos de simulação nu- 

mérica de colectores so- 

lares planos ... «e 444 
E. Maldonado 
E. Oliveira Fernandes 


11-977 65 


659.8 — PROTECÇÃO , IMPERMEABILIZAÇÃO E 
ISOLAMENTO DAS CONSTRUÇÕES 


C. D. U. 699.81 


Acção dos incêndios sobre 

& Construção ... cc ves 445-446  12-1-977/78 231 
Alfredo Manuel Filorenço Tovar de Lemos 
Illefonso Cabrita Neves 
Joaquim da Conceição Valente 


72 — ARQUITECTURA 


C. D. U. 72.011.1 


Necessidades humanas e 
exigências funcionais 
da habitação ... ... 

Ruy José Gomes 


445-446  12-1-977,/78 191 


C. D. U. 72.011.1:69.002 


Tecnologias da construção 
de habitação ... .. .. 448 
José A. Teixeira Trigo 


3-978 463 


a técnica moderna emprega 


BETÕES REFRACTARIOS 


à base dos cimentos 
FONDU LAFARGE 
SECAR 250 


porque são 
práticos 
eficientes 
econômicos 


REFRACTÁRIOS ATÉ 1800º C 
REFRACTÁRIOS ISOLANTES ATÉ 1700º C 
Os nossos serviços técnicos, especializados, estão à vossa inteira 


disposição para estudar, sem qualquer encargo da vossa parte, a 
adaptação dos BETÕES REFRACTÁRIOS a todos os vossos problemas 


AGUIAR & MELLO LL.” 
P. do Municipio, 13-1.º — LISBOA — Tel. 321151/2 


TÉCNICA XIII 


À nossa vocação é exportar 


Caixa de escape para uma central nuclear, 
destinada à Electricité de France 


EQUIPAMENTOS HIDROMECÂNICOS 


EQUIPAMENTOS ELECTROMECÂNICOS 
PARA CENTRAIS HIDROELÉCTRICAS 

E TÉRMICAS 

(CLÁSSICAS E NUCLEARES) 


EQUIPAMENTOS PARA A INDUSTRIA 
QUÍMICA E PETROLÍFERA 


MATERIAL CIRCULANTE DE CAMINHO 


DE FERRO Países para onde temos exportado : 
EQUIPAMENTO OFFSHORE PARA AFRICA DO SUL - ANGOLA - ARGENTINA - BRASIL - COSTA 
E e PE ARA BAEIO RICA - FRANÇA - GUATEMALA - IRÃO - IRAQUE - ISLÂNDIA 
DE PETROLEO KUWEIT - MALAWI - MARROCOS - MÉXICO - MOÇAMBIQUE 
CONSTRUÇÕES. MODULADAS, NICARÁGUA - PAQUISTÃO - RODÉSIA - SUDÃO - TURQUIA 
PONTES E ESTRUTURAS UGANDA - USA. - VENEZUELA - ZAIRE 


SOREFAME 


SOCIEDADES REUNIDAS DE FABRICAÇÕES METÁLICAS, S.A.R.L. 
Amadora — Portugal 


Edições da Técnica 


TABELA DE PREÇOS 


TABELAS PARA O CALCULO DO BETÃO ARMADO 


Brazão Farinha, 6.' edição — 1970 


TOPOGRAFIA GERAL 


C. Xerez— 1.º vol. 3.º edicão — 1978 


C. Xerez— 2.º vol. 2,º edição — 1966 


TABELAS TECNICAS 


Vicente Ferreira e Brazão Farinha, 8.º edição 


1977 


MANUAL DE HIDRÁULICA GERAL 


A. Lencastre, 2.º edição — 1969 


DICIONÁRIO DE UNIDADES E TABELAS 
DE CONVERSÃO 


Vasco Costa e Osvaldo Francês, 1.º edição — 


1959 


NUMEROS E VARIÁVEIS COMPLEXAS 
A. H. Simões de Abreu, 1.º edicão 1977 


4508C0 


350800 


250500 


500800 


450800 


75800 


312550 


(*) LIVROS A PUBLICAR BREVEMENTE 


Desconto de 10% aos assinantes 


GUIA DE ANALISE QUÍMICA DAS ÁGUAS 


A. Herculano de Carvalho, 1.º edição — 1961 


CALCULO DE PÓRTICOS 


G. Kani — Trad. da 7.º edição alemã — 1962 


MANUAL DE ESTRUTURAS 
Brazão Farinha — 1,º VOLUME 
Encadernado 


Plastificado 


2º e 3.º Volumes 


BETÃO PREESFORÇADO 
(A cargo dos participantes na Semana de Be- 


tão Preesforçado — |. S, T.— 1975) 


150500 


100800 


625$00 


587850 
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425400 
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